Collisional properties of ultracold K-Rb mixtures 
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We determine the inter-species s-wave triplet scattering length as for all K-Rb isotopic mixtures 
by measuring the elastic cross-section for collisions between 41 K and 87 Rb in different temperature 
regimes. The positive value £13 = +1631*2 a ° ensures the stability of binary 41 K- 87 Rb Bose-Einstein 
condensates. For the fermion-boson mixture 40 K- 87 Rb we obtain a large and negative scattering 
length which implies an efficient sympathetic cooling of the fermionic species down to the degenerate 
regime. 
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Ultracold mixtures of alkali atoms have recently been 
the subject of intensive experimental research, mainly 
aimed at the production of novel degenerate quantum 
systems. Sympathetic cooling 1] in a mixture represents 
a way to reach the degenerate regime for species which 
cannot be efficiently cooled by direct evaporation, as in 
the case of 6 Li and 41 K In addition, mixing differ- 
ent isotopes or atomic species allows to create ultracold 
Bose- Fermi systems 0, Q , and paves the way to the cre- 
ation of binary Bose-Einstein condensates with different 
atomic species and of ultracold heteronuclcar molecules. 

In an ultracold mixture the interaction between the 
components is described in terms of a single key parame- 
ter, the inter-species s-wave scattering length, which for 
example determines the efficiency of the sympathetic- 
cooling procedure. In addition, knowledge of the scatter- 
ing length is important in order to predict macroscopic 
properties of these mixed systems, such as the stability 
and phase separation of binary degenerate systems, or 
microscopic properties, such as the occurrence of Fes- 
hbach resonances in the inter-species collisions. While 
collisions either between alkali of the same species or be- 
tween different isotopes have often been well character- 
ized [jj, collisions between different alkali species repre- 
sent a still relatively unexplored field. 

In this context, K-Rb mixtures are particularly inter- 
esting, since they are promising candidates in the gen- 
eral field of Bose-Fermi and Bose-Bose degenerate gases. 
Light-assisted inelastic collisions between various K-Rb 
isotopes have been investigated in magneto-optical traps 
(MOTs) 00. Ultracold elastic collisions have been in- 
stead studied in a magnetic trap, but only the magnitude 
of the 41 K- 87 Rb s-wave triplet scattering length has been 
estimated 3]. In this Letter we report the results of a 
thorough investigation of the collisional properties of a 
41 K- 87 Rb mixture. By studying the interspecies elastic 
collisions in a broad temperature range, we determine 
the sign of the scattering length and tighten the bounds 
on its magnitude. Thanks to this result, we are able to 
derive by mass scaling the triplet scattering lengths for 
all the other combinations of K and Rb stable isotopes 



and to discuss the possibility of producing mixed K-Rb 
degenerate systems. We also determine a region of confi- 
dence for the singlet s-wave scattering length of the 41 K- 
87 Rb mixture from inelastic decay rates of a mixed-spin 
sample. 

The ultracold 41 K- 87 Rb mixture is produced using the 
apparatus described in In brief, about 10 7 K atoms 
at 300 [iK and 5 x10 s Rb atoms at 100 /xK are loaded in 
a QUIC magnetostatic trap using a double-magneto- 
optical trap scheme. Both species are prepared in their 
doubly-polarized spin state \F = 2, Mp = 2) thanks to an 
optical pumping stage. A gas prepared in such state has 
the advantage of being immune to spin-exchange losses 
and is typically stable. Evaporative cooling is then per- 
formed selectively on the Rb sample using a microwave 
knife tuned to the hyperfme transition at 6.8 GHz (the 
corresponding transition for K is at 254 MHz), while 
the K sample is sympathetically cooled through elastic 
K-Rb collisions. The number of atoms and tempera- 
ture of both samples is determined simultaneously at 
the end of each experimental run by absorption imag- 
ing, using two delayed resonant light pulses after release 
from the magnetic trap. Both atomic species experience 
the same trapping potential with cylindrical symmetry. 
The axial and radial harmonic frequencies for Rb are 
ojax = 2ttx 16.3 s" 1 and ui ra d — 2ttx 197 s _1 respectively, 
while those for K are by a factor (M Rb /M K )^ 2 = 1.45 
larger. 

We have investigated the dependence of the K-Rb elas- 
tic cross-section on temperature, by performing a series 
of rethermalization measurements. We drive the system 
out of thermal equilibrium and observe the subsequent 
thermalization mediated by elastic collisions. The differ- 
ent trap frequencies for the two species allow us to selec- 
tively heat Rb by superimposing a 10% modulation on 
the trapping potential at twice the Rb radial frequency 
for typically 100 ms. As an example we show in Fig. Q] 
the measured increase of temperature of the K sample at 
T « 1.6 /iK that follows heating of Rb, together with an 
exponential fit to the data. The time constant for rether- 
malization t is determined from the fit. We have checked 
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that the aspect ratio of the K sample was constant dur- 
ing the whole rethermalization which is evidence for a 
common temperature of the spatial degrees of freedom. 

We plot in Fig. [21 the measured elastic-collision rate 
l/(rfi) as a function of temperature in the range (1.6 — 
45) fiK. Here n = (-fi— + T^) J riKnRbd 3 x is the effective 
density of K and Rb atoms. A general model to describe 
the thermalization between atoms of different masses can 
be obtained by extending the one of Ref. |h| t o energy- 
dependent cross sections and to partial waves higher than 
s-wave. We have checked that at our temperatures only 
s- and p-wave collisions are relevant. The thermalization 
time per unit atomic density mediated by s- and p- wave 
collisions is 



1 



2a s k B T 



(a,vE) + —(a P vE) 



(1) 



where a s ~ 2.7 is the average number of s-wave collisions 
necessary for rethermalization, v is the relative velocity of 
a colliding pair and E is the relative collision energy. We 
obtain the ratio ^ = f upon integration over all collision 
directions. The factor £ = 4/x/M, with M and [i the 
total and the reduced mass respectively, gives a reduction 
£ ~ 0.87 of the thermalization efficiency with respect to 
the case of equal masses. Finally, the averages (■) are 
performed on a classical Maxwell-Boltzmann distribution 
of relative velocities. 

The partial cross sections are obtained from a numeri- 
cal solution of the scattering equations using standard 
propagation algorithms. The molecular Hamiltonian 
for collisions of atoms with hyperfine structure includes 
the kinetic energy of the relative motion, the hyperfine 
atomic energy and the adiabatic Born-Oppenheimer 
and 3 S+ symmetry potentials. The latter have been con- 



structed by smoothly matching the short-range ab ini- 
tio potential of |ll| onto a long-range analytic potential 
V = V d ± Vex/2. Here V d = -C 6 /R 6 ~ C 8 /R 8 - C w /R 10 
is the dispersion potential and V ex is the exchange poten- 
tial. An accurate van-der-Waals coefficient C§ is taken 
from Ref. C*s and Cio from ^| an d the analytic 

form of V ex from [lj|- The molecular potential is made 
flexible by adding a short-range correction to the adia- 
batic potentials (see Ref. 0). This procedure allows 
us to tune singlet and triplet scattering lengths to agree 
with the data. 

The experimental data are fitted using a min-% 2 pro- 
cedure. Since collisions between atoms prepared in a 
doubly-polarized spin state are single-channel collisions 
involving only the 3 S+ symmetry potential, the scatter- 
ing length d3 is the only free parameter in the model |16| . 
Our main results are shown in Fig. [21 The best agreement 
is obtained for a positive scattering length 03 = 163 do, 
while the best-fit curve for negative 0,3 fails to fit the 
experiment at high temperature (solid lines in Fig. 
Actually, for 03 < the cross-section drops with energy 
from its threshold value a s — Aira^ at lower collision en- 
ergies than it does for 0,3 > 0. Moreover, the curve for 
as = 163 ao has a significant contribution from a broad 
p-wave shape-resonance near the top of the centrifugal 
barrier which further increases the rate at high T. This 
resonance rapidly shifts at lower energies for decreasing 
0,3 (dashed line in Fig. This circumstance and the 
absence of observed resonant features sets a tight lower 
bound 03 w 150 ao on the confidence interval for 0,3, while 
the upper bound is looser. Actually, scattering lengths 

~ 200 ao having only a minor p-wave contribution still 
agree well with the data. 

In order to confirm our determination of the magni- 
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FIG. 1: Evolution of the mean temperatures of the K sample 
after 100 ms of selective heating of Rb. Each datapoint is the 
average of five measurements, and the error bar displays the 
corresponding statistical fluctuation. The continuous line is 
the best fit to an exponential function; the time-constant is 
inversely proportional to the K-Rb elastic collisional rate. 



FIG. 2: Dependence on temperature of the elastic collision 
rate l/(rn). The points are the experimental data. The solid 
lines are the best fit to the model described in the text for 
positive and negative triplet scattering lengths, corresponding 
to C13 = 163 ao and 03 = —209 ao, respectively. The dashed 
line corresponds to 03 = 140 ao. 
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tude of 03 we have performed an independent measure- 
ment of the collisional cross-section, studying the damp- 
ing of dipole oscillations of a sample composed by 8x 10 4 
K atoms and 1.5xl0 5 Rb atoms, at T = 1.7/iK. Dipolar 
modes of the two atomic clouds are excited by displacing 
the potential minimum along the weak horizontal trap 
axis for about 30 ms. As shown in Fig. [3] the center 
of mass oscillations get rapidly coupled and damped by 
elastic collisions. We fit the experimental data to the so- 
lutions of the coupled equations of motion for the clouds 
centers of mass 
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where the labels 1 and 2 denote K and Rb respectively, 
M = Mi + M 2 , N = Ni + N 2 , T = na s v is the s-wave 
collision rate, v is the rms relative velocity, and n is the 
atoms density defined above. The low temperature at 
which the experiment is performed ensures that s-wave 
collisions dominate. Our system appears to be in the 
collisionless regime ( V << u;i,<j 2 ), where the position 
of each center of mass is the sum of two damped sinu- 
soids at the bare frequencies of K and Rb. Our best fit, 
shown in Fig. [21 well reproduces the experimental ob- 
servation. From the fitted collision rate T = 5.2(7) s _1 
we determine the magnitude of triplet scattering length 
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I70I35 a o> where the uncertainty is dominated by 
that on the atom number. This value for I03I is in good 
agreement with the one previously derived from thermal- 
ization measurements. 

We determine the singlet K-Rb scattering length a±, by 
measuring the trap loss for collision between K atoms in 
the hyperfine state |2, 2} and Rb atoms in |2, 1). Our sam- 



ple is initially composed by a small number of K atoms 
(3xl0 4 ) and a larger number of Rb atoms (1.2xl0 5 ), 
both in |2,2), at a temperature T = 1.8 /iK. We then 
transfer a small fraction (20%) of both species in |2, 1) us- 
ing a radio-frequency sweep and we study the subsequent 
decay of K, which takes place on a timescale much shorter 
than the lifetime of the doubly spin-polarized mixture. In 
these experimental conditions K atoms in |2, 1} and Rb 
atoms in |2,2) do not contribute significantly to loss of 
K atoms from the trap. The time evolution of K-atoms 
number is given to a good approximation by the solution 
of the single rate equation h K {t) = — Gi 



K 

2,2) 



(<)^,)(t) 
Here 



with the constraint N^ 2 ^(t) — NS^^t) = const. 
n is the spatial density, N the atom number and G is 
the rate constant for inelastic collisions between K in 
|2,2) and Rb in |2, 1} \r%. Our collision model agrees 
well with the experimental data for an inelastic collision 
rate G=1.8(9)xl0" 11 cm 3 s _1 . A comparison between 
the measured inelastic rate and the numerical calcula- 
tion (shown in Fig. 0J) allows us to determine the region 
of confidence for the singlet scattering length a\, given 
our interval of 0,3 values. This comparison singles out a 
region of positive singlet scattering length, a\ > 30 ao, 
excluding a broad region of strong suppression of inelastic 
processes centered at a\ ~ 03. In alternative, the large 
negative singlet scattering length a\ < —210 ao is also 
compatible with the data. This result represents a first 
characterization of the singlet K-Rb interaction, and can 
provide a useful guidance for further collisional studies, 
in particular for the investigation of Feshbach resonances 
in selected spin-states. 

We have also determined the triplet s-wave scatter- 
ing lengths for all the other pairs of K and Rb isotopes. 
It must be remarked that our analysis for 41 K- 87 Rb is 
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FIG. 3: Dipole oscillations of the K (a) and Rb (b) clouds 
along the weak trap axis. The solid lines are a best fit to 
the model presented in the text. The faster damping and the 
beatings observed in the K motion are due to the smaller mass 
of the K sample. 



FIG. 4: Contour plot of the calculated inelastic collisional 
rate G in units of 1CP 11 cm 3 s _1 vs the triplet scattering 
length 03 and the reduced singlet scattering length x\ = 
2 /it arctan(ai/a sc ), where a sc = 72 ao. The shaded regions 
correspond to the rate determined in the experiment. 
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TABLE I: Triplet s-wave scattering lengths (in a.u.) for colli- 
sions between K and Rb isotopes computed using the nominal 
triplet potential. 

39 K 40 K 41 K 



Rb 58±l 4 -38+H 3291 1 ™ 
Rb 3ltg 6 -261^12 163±g 



weakly dependent on the number Nb of bound states sup- 
ported by the ABO potentials. In contrast, mass scaling 
from one isotope to the others depends more strongly on 
Nf,. Tab. H] summarizes our results for the number of 
bound states N{> = 32 supported by the nominal triplet 
potential. Calculations that allow for a realistic (±2) un- 
certainty in Nf, do not qualitatively change the results of 
Tab. [I] With this proviso, we discuss below the conse- 
quences of these results. The repulsive character of the 
41 K- 87 Rb triplet interaction indicates that it is possible 
to form a stable binary Bose-Einstein condensate. In- 
deed, in case of a negative inter-species scattering length 
much larger in magnitude than the single-species ones 
[l8| . the mean- field interaction would lead to a collapse 
of the condensates 0, . Of great interest is also the 
40 K- 87 Rb p a j rj f or which we g nc i a negative and probably 
large inter-species scattering length. Therefore, at least 
in the absence of unfortunately placed zeros in the elas- 
tic cross section, sympathetic cooling of fcrmionic K with 
87 Rb can be expected to be highly efficient. Moreover, 
the attractive character of the interaction should prevent 



phase separation of the two components once the degen- 
erate regime has been reached |23j . thus insuring an effi- 
cient thermalization even in this regime. This represent a 
qualitatively different situation with respect to the other 
known Bose-Fermi mixture, 6 Li- 7 Li, for which the inter- 
species scattering length is positive FTil. We note that 
among the pairs containing 85 Rb, which however does 
not seem to be an easy partner for sympathetic cooling 
[25|], only the 41 K- 85 Rb pair shows a significantly large 
interaction. Finally, the mixture 39 K- 87 Rb is not partic- 
ulary appealing, since the interspecies scattering length 
is expected to be relatively small and, moreover, the for- 
mation of a stable Bose-Einstein condensate of 39 K with 
a large number of atoms is prevented by the negative sign 
of its own scattering length . 

In conclusion we have drawn a picture of the triplet 
interaction between all the K-Rb isotopes pairs. Our 
study paves the way to the production of novel mix- 
tures of degenerate atomic gases, such as Bose-Bose sys- 
tems with large repulsive interaction and Bose-Fermi sys- 
tems with attractive interactions. Moreover, it provides 
important information for predicting the occurrence of 
magnetically-tunable Feshbach resonances. These reso- 
nances would represent a tool to control the inter-species 
interactions and to drive the production of ultracold het- 
eronuclear molecules. 
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